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Abstract 



Using an interface displacement model we calculate the shapes of thin 
liquidlike films adsorbed on flat substrates containing a chemical stripe. We 
determine the entire phase diagram of morphological phase transitions in these 
films as function of temperature, undersaturation, and stripe width. 
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I. INTRODUCTION 



At present various experimental techniques allow one to provide substrates with well- 
defined, microscopically small geometrical or chemical structures or combinations thereof 
(see, e.g., Ref. These structured substrates can serve as devices to guide and process 
tiny amounts of liquids. Such "microfiuidics" systems are promising for future appli- 
cations, e.g., in chemistry and biology [^3- Although the performance of these systems 
ultimately depends on their dynamic properties, a thorough understanding of the corre- 
sponding structural properties in thermal equilibrium is an important prerequisite [^. In 
this spirit as a paradigmatic case we determine the equilibrium structures of liquidlike wet- 
ting films forming at the interface between a bulk vapor phase and a fiat solid substrate 
containing a single chemical stripe (Fig. p. Our aim is to explore the entire phase dia- 
gram for the emerging lateral fiuid structures as function of temperature, undersaturation 
or pressure, and stripe width. 

In Ref. it has been shown that a thin liquidlike film adsorbed on such a substrate 
(Fig. |l|) can undergo an interesting morphological phase transition (Figs. |^ and ^) of the 
shape l{x) of the liquid- vapor interface, depending on the structures of the effective interface 
potentials A-|-(/) characterizing the corresponding homogeneous substrates and "— " 
(see, e.g.. Fig. ^). In accordance with Ref. P we study a substrate designed such that the 
effective interface potential A+(/) of the stripe part exhibits two competing local minima 
("I" near and "III" further away from the substrate surface) whereas the asymptotic film 
thickness Z_ far from the stripe is given by the single global minimum ( "II" ) of the interface 
potential A_(Z) of the embedding substrate (see Fig. H). In Ref. the phase diagram for 
this morphological phase transition has been determined, at constant temperature T, as 
function of the stripe width a and of the undersaturation A/i = /^o ~ > where fioiT) is 
the chemical potential at liquid-vapor coexistence. The purpose of the present study is to 
determine how this phase diagram for the morphological phase transition within the a-A^ 
plane evolves as function of temperature. This additional information is important in order 
to be able to put the morphological phase transitions into the context of T-Afj, wetting phase 
diagrams, including the prewetting line, on homogeneous substrates 

We emphasize that experimental techniques such as refiection interference contrast mi- 
croscopy (RICM) and atomic force microscopy (AFM) in tapping mode are capable of 
imaging the liquidlike structures under consideration down to the nm scale. Moreover, due 
to their relatively small spatial extensions the structures discussed here are well accessible 
by Monte Carlo simulations as another means for testing the theoretical predictions. For 
these envisaged studies the knowledge of the full topology of the phase diagram as function 
of T, A/i, and a is very important. 



We use the following simple interface displacement model within which the equilibrium 
interface profile l{x) minimizes the functional 
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where A = L^Ly is the area of the flat substrate surface located at 2; = 0. The flrst part takes 
into account the effective interaction between the hquid-vapor interface and the substrate 
via the laterally varying effective interface potential A(x, /) = Aflhl + ijj{x,1). The second 
part in Eq. (P is the leading-order term of a gradient expansion of the surface free energy 
associated with the deviation of the liquid- vapor interface profile from its flat conflguration; 
(7ig is the surface tension. AQb = '^fJ'iPi — Pg) + is the difference of the bulk 

free energy densities of the liquid and the vapor phase with number densities pi and pg, 
respectively. Alternatively, in terms of the bulk pressure p, AVLh is given approximately by 
the Gibbs-Thomson expression ksTpi \ia{psat/p) where Psat is the saturated vapor pressure. 
At liquid-vapor coexistence A{x,l) is given by uj{x,l) which is of the form Yli>2^ii^)/^^ ■ 
For reasons of simplicity we assume that the lateral variation of uj{x, I) is steplike: uj{x, I) = 
6(|x| — a/2)ci;_(/) -|- 0(a/2 — |a;|)a;+(/) with the effective interface potential (c<j_) of a 
homogeneous substrate composed of ("— ") particles. In Refs. [|lOl and [jlT[ we have 
shown that the functional in Eq. (|I|) - despite its simplicity - allows one to determine 
reliably the morphology of liquidlike films on structured substrates, and we have presented 
a derivation of this square-gradient functional from a microscopic density functional theory. 
We also found (see Fig. 11 in Ref. |T^) that the aforementioned steplike variation of uj{x, I) 
as opposed to its actual smooth lateral variation represents a rather reliable approximation. 
The form of uj{x,l) follows from considering as the interaction potential between the fiuid 
particles a Lennard- Jones potential with a depth — e and an effective particle diameter a. 
Moreover, also the interaction between the fluid and the substrate particles is modeled 
by a Lennard- Jones potential (compare Ref. HT^])- Since the substrate is translationally 
invariant in the y direction the effective interface potential depends only on x, so that l{x) 
is also translationally invariant in the y direction. The interaction potential parameters 
of the stripe part ("+") are the same as for the "+" substrate described in Fig. 11 in 



Ref. |]I0[, and for the embedding "— " substrate part we use the parameters for the 
substrate part in Fig. 7 in Ref. ||T^ multiplied by 0.9102. The temperature dependence 
of the effective interface potentials A-i-(Z) is determined by systematically keeping track of 
the temperature dependence of the bulk liquid and vapor densities pi and Pg which enter 
into A-|-(/) as discussed in Sec. II in Ref. |]10[- For the case studied here the temperature 
dependence of uj+{l) for the stripe part is such that a homogeneous substrate exhibits 
a first-order wetting transition at T* = ksTyj/e ~ 1.102 and prewetting transitions along 
the prewetting line AppreiT), T > and Ap > 0, that extends into the vapor phase 
region of the phase diagram (thick dashed line denoted as "P" in Fig. |^). ui-{l) for the 
embedding substrate is also temperature dependent such that the asymptotic film thickness 
grows with increasing T and shrinks with increasing Ap. As it turns out, however, the 
exact temperature dependence and the magnitude of /_ are not important for the qualitative 
features of the phase diagram for the morphological transitions. 

The systems under consideration here differ from those studied in Ref. |]T2[ in two im- 
portant aspects. First, here the liquidlike films are so thin that they are completely under 
the infiuence of the substrate potential which determines, inter alia, the effective interface 
potential. Secondly, we consider a grand canonical ensemble without a volume constraint 
for the liquid phase so that the liquidlike films are not subject to instabilities along the y 
direction as observed in, e.g., Ref. [|1^ . 



Using a numerical relaxation technique |TB[ we solve the two-point boundary value prob- 
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lem for the Euler-Lagrange equation 

with l{x ±00) = /_ (2) 



d^'^ l=l(x) 

which follows from Eq. (|1]) by functional differentiation with respect to l{x). Equation (0) 
exhibits the structure of a one- dimensional classical mechanical equation of motion in a 
time-dependent external potential. With the stripe width a and the temperature T fixed, 
for a wide range of undersaturations there are two solutions of Eq. which correspond 
to local minima of the functional VLs- One of the solutions is closely bound to the stripe 
whereas the other solution is further away or even repelled from the stripe as described in 
Ref. 0. Figures show pertinent examples for the effective potentials K± of the two 
substrate materials and the two appertaining solutions of Eq. (H) at three different points 
of the phase diagram: Figs. ||, ||, ^ and |^ correspond to the points (T, A/i) indicated by 
"a", "b", "c", and "d", respectively, in the phase diagram shown in Figs. ^ and 0. The 
equilibrium solution l{x) for the profile corresponds to the global minimum of ^5, and the 
other solution is metastable. At a certain value A/i^, for which Vlg has the same value for 
both solutions, a phase transition from one interfacial configuration to the other takes place. 
This transition is first order because the derivatives of VLs{l{x)) with respect to A/i are 
discontinuous. 

As already shown in Ref. in Eq. (|I|) the two contributions to Vtg compete with each 
other in minimizing the whole functional. Depending on the special choices for a, T, and 
A/i the cost in free energy for increasing the liquid-vapor interface area is overcompensated 
by the gain in free energy which follows from occupying the deeper minimum (I) of A+ near 
the wall, so that the bound configuration has the lower free energy. In the opposite case 
this gain in free energy is too small to compensate the cost in free energy from the increased 
area of the liquid-vapor interface, leading to a less bound or even repelled interface profile. 



III. DISCUSSION OF THE PHASE DIAGRAM 

In Fig. the lines A/it(T; a) = ftoiT) — UtiT] a) of phase coexistence between the bound 
and the repelled configuration for different stripe widths a are presented. These lines are 
vertical cross-sections of the full phase diagram shown in Fig. |^. For a given stripe width a the 
bound (repelled) solution is stable for —A/i < — A/i((T; a) (—A/i > — A/it(T; a)), i.e., below 
(above) the corresponding line of phase coexistence. The triangles indicate the intersections 
To(a) between the coexistence lines A/if (T; a) and the bulk liquid-vapor coexistence line 
A/i = 0. At To (a) the morphological phase transition occurs at liquid- vapor coexistence, 
forming the line denoted as a in Fig. |^. Upon increasing the stripe width a the lines of 
coexistence are shifted towards the prewetting line "P" without touching or crossing it for 
any finite value of a. (This means that for large a the lines of coexistence do not end and 
reappear at the prewetting line such as, e.g., the prefilling transition lines in a wedge-shaped 
groove for large opening angles of the wedge |1^.) For constant undersaturation A/i we find 



Tpre(A/i) — Tf(A/i;a — >■ cxo) ~ a~ for the difference between the transition temperatures 
Tt(A/i;a) for the morphological and the prewetting transition. In the limit a ^ 00 the 
morphological phase transition tallies with the prewetting transition on the homogeneous 
substrate. 
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On the mean-field level considered hitherto each coexistence curve Afit{T; a) is a line of 
first-order transitions ending in a critical point (Tc(a), A/ic(ci)) (denoted as full diamonds in 
Fig. 0and forming the line denoted as 7 in Fig. |) such that at T > Tc(a) or —Afj, < — A/ic(o) 
for the given stripe width a there is no morphological transition but a smooth variation from 
the bound to the repelled solution. As demonstrated by Fig. ^ and its inset the positions of 
the critical points exhibit a nontrivial and nonmonotonous dependence on the stripe width 
a; there is no simple criterion for the effective interface potential A± which allows one to 
predict the corresponding line of critical points (see the line 7 in Fig. 0) as function of a. 
As already pointed out in detail in Ref. |^ interface fluctuations along the y direction in 
this effectively one-dimensional stripe configuration actually smear out the sharp first-order 
morphological phase transitions and thus eliminate the critical points [l^. In Fig. |^ for 
a = 5a the width of this fluctuation-induced smooth transition region is indicated by the 
thin dashed-dotted lines. For a > 10a these fluctuation effects are already negligibly small, 
apart from the close vicinity of the critical points which are still eliminated for any finite a. 
However, in the limit a — ^ 00 the coexistence lines Afit{T; a) merge with the prewetting line 
"P" (see Fig. |^) associated with the homogeneous "+" substrate, which does have a genuine 
critical point "C" beyond mean field theory. 

The three thermodynamic states denoted as "a" , "b" , and "c" in Figs. ^ and ^ lead to the 
interface configurations shown in Figs. ^ |^, and |^, respectively. The corresponding effective 
interface potentials A-|-(/) differ with respect to the number and the relative positions of their 
local minima. For temperatures T < Tpre and small undersaturations A/i, A+(Z) exhibits 
two local minima (I and III); an example for this case is shown in Fig. ^ appertaining to 
the thermodynamic state denoted as "a" in Fig. 0. This case is analogous to that described 
in Ref. |0. At the thin dotted line A/ij(T) denoted as "i" in Fig. |^ the minimum III 
of A_|_ far from the wall and the maximum of A+ merge, forming a saddle point; A/ij(T) 
ends at the prewetting critical point "C". For — A/x < — A/ij(T) there is only one local 
minimum (I) of A_|_. An example for this latter situation is shown in Fig. ^ corresponding 
to the thermodynamic state denoted as "b" in Fig. 0. Surprisingly the coexistence lines 
Afit{T] a) for the morphological phase transition extend below this line "i". This shows that 
the presence of an energy barrier and of a second local minimum of A_|_ is not a necessary 
condition for the occurrence of the morphological phase transition described here and in 
Ref. 0. The thin dotted line AumiT) denoted as "iii" is the line at which, analogously, 
the minimum I near the wall and the maximum of A+ merge, forming a saddle point. Both 
lines "i" and "iii" meet and end at the critical point "C" of the prewetting transition of a 
homogeneous "+" substrate where both minima I and III and the maximum of A_|_ merge, 
leaving a single minimum. 

At the thin short-dashed line A^uiT) denoted as "ii" in Fig. ^ the positions of the 
second minimum III of A+ and of the minimum II of A_ coincide. Along this line one of the 
solutions of Eq. (^ is completely flat: due to dA±/dl (/ = /_) = Eq. (0) leads to the trivial 
solution l{x) = /_ with (Pl{x)/dx'^ = 0. An example for this case is shown in Fig. ^ which 
corresponds to the thermodynamic state denoted as "c" in Fig. |^. The line "ii" separates 
the region where the repelled solution is bent towards the stripe (— A/i < —AnuiT)) from 
that where it is bent away from the stripe (— A/x > —AfiaiT)). We note that the fact that 
one of the solutions is completely flat is an artifact of our simplifying model assumption 
that u!{x, I) varies steplike as function of x. In a more realistic model with a smooth lateral 
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variation of the effective interface potential (as used in Ref. 0) tlie corresponding solution 
would be almost fiat with a small curvature. In the system considered here the asymptotic 
film thickness /_, i.e., the position of minimum II, is always larger than the position of 
minimum I. However, for low temperatures there is a line at which the minimum II of A_ 
and the maximum of A+ coincide. This occurs along the thin long-dashed line A/ij^,(T) 
denoted as "iv" in Fig. This line is the continuation of line "ii" for lower temperatures T 
starting at the point where the line "ii" meets the line "i" tangentially. Along this line, too, 
the repelled solution of Eq. is constant, i.e., l{x) = Z_. As the line "ii", also line "iv" 
separates two regions with different curvature behavior of the repelled solution. 

For temperatures T > Tp^e, i.e., on the right side of the prewetting line P, the minimum 
III of A+ is deeper than the minimum I. The difference between the asymptotic film thickness 
/_ corresponding to the minimum II and the deeper minimum III is also smaller than that 
between /_ and the local minimum I near the wall. Therefore both contributions to Qs 
are larger in the case that the interface follows minimum I as compared to the case that it 
follows minimum III so that for T > Tp^e the bound solution is always metastable. This 
situation is shown in Fig. |^ which corresponds to the thermodynamic state denoted as "d" 
in Fig. 1^. If T > TmlAfi) A+ exhibits only one local minimum (III) and therefore there is 
only one solution of Eq. (^, namely the repelled one. 

IV. CONCLUSION 

Within an interface displacement model based on a microscopic density functional theory 
we have calculated the entire phase diagram (Figs. |^ and ^ of morphological phase transi- 
tions of wetting films on a substrate with a chemical stripe (Fig. |l]). This phase diagram and 
the corresponding equilibrium interface profiles shown in Figs. ^-|^ elucidate the dependence 
of the morphological phase transitions on the thermodynamic variables temperature and un- 
dersaturation (or, equivalently, pressure) as well as on the stripe width. The morphological 
phase transitions have been put into the context of the prewetting transitions occuring on 
a homogeneous substrate formed by particles of the stripe material. 
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FIGURES 
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FIG. 1. Schematic cross-section of the morphology of a hquidhke film covering a planar 
substrate which contains a slab of different material. In top view the substrate forms a stripe 
composed of "+" particles extending from a; = —a/2 to a; = a/2. The substrate is translationally 
invariant in the y direction, with its surface located at z = 0. Z_ = l{\x\ —i- oo) is the equilibrium 
film thickness corresponding to a homogeneous substrate composed of "— " particles. The stripe 
(embedding) part exerts an effective interface potential A_|_ (A_) on the liquid- vapor interface at 
the vertical distance z = l{x). In the case shown here the "— " substrate prefers a thicker wetting 
film than the substrate. 
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FIG. 2. T-A/i phase diagram for the morphological transitions of a liquidlike layer on a 
chemical stripe (compare Fig. |l|); T* = ksT/e and /i* = /i/e. A/i = is the bulk liquid- vapor 
coexistence line. For the present model the bulk critical temperature is T* ~ 1.412; the estimated 
triple point temperature is ?a 0.8. The homogeneous substrate composed of "+" particles, i.e., 
corresponding to the stripe, supports a first-order wetting transition at coexistence denoted by 
"W" (0). The corresponding prewetting line "P" (thick dashed line) is attached tangentially and 
ends at a prewetting critical point "C" (□). The effective interface potentials A± and, as examples, 
the equilibrium interface profiles at the thermodynamic states • denoted as "a" , "b" , "c" , and "d" 
are shown in Figs. |^-|^ for a/ a 27.5, 6.4, 26, and a/ a = 50, respectively. For "a", "b", and "c" 
these values for the width a are those for which in these thermodynamic states the morphological 
phase transition takes place. The coexistence lines A^j(T;a) (full lines) of the bound and the 
repelled solutions for fixed stripe width a are accompanied by the appertaining value a/ a. For 
— A|U > —Ant the repelled solution and for — A^ < —A/it the bound solution is the stable one. 
The lines of coexistence intersect the axis Afi = at the points To(a) marked as ▲. Afit{T;a) 
joins the bulk liquid-vapor coexistence curve tangentially with Afit{T —>■ ?b(a); a) ~ (T — To(a))^; 
the numerical data indicate that 6 is larger than the corresponding exponent 5 = 3/2 for the 
prewetting line. Due to this tangential behavior the intersections ro(a) occur at lower temperatures 
than expected visually on the present scale. Therefore we have separately labeled the intersections 
To (a) by the corresponding values for a/ a. All lines of coexistence end at a critical point denoted 
by The dashed-dotted lines accompanying the line for a = 5a show the region within which 
the fluctuation-induced rounding of the first-order phase transition takes place; for a > lOo" this 
rounding effect is not visible on the present scale. The inset magnifies the region around the 
prewetting critical point "C" . The meaning of the lines "i" , "ii" , "iii" , and "iv" is discussed in the 
main text. 
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FIG. 3. The thermodynamic states at which the morphological phase transitions take place are 
located on the surface "S" . The lines of coexistence shown as full lines in Fig. ^ are slices through 
this surface "S" at constant values of 1/a. "S" is a surface of first-order transitions and separates the 
region where the repelled solution is stable (above "S" ) from the region where the bound solution 
is stable (below "S"). It is bounded by the prewetting line (dashed line "P") which lies in the plane 
1/a = 0, by the line 7 (dotted line) of critical points of the morphological phase transitions, and 
by the line a which lies in the plane A/x = indicating the loci (To(a), 1/a) of the morphological 
phase transitions at bulk liquid-vapor coexistence. "W" and "C" indicate the first-order wetting 
transition on a homogeneous substrate (0) and the prewetting critical point (□), respectively; 
both lie in the plane 1/a = 0. Indicated by the thin line in the front, "S" extends out to larger 
values of 1/a where, however, it is increasingly smeared out by the fluctuation-induced rounding 
of the morphological phase transition (compare Fig. |2[). Due to Tpre(A;u) — Tt(A/i; a — > cxd) ~ a~^, 
in terms of 1/a "S" is approximately a ruled surface: one of the two principal curvatures is very 
small, vanishing exactly in the limit 1/a 0. The thermodynamic states denoted by "a", "b" 
and "c" (•, compare Fig. ^) lie on the surface "S" whereas the thermodynamic state "d" is located 
above "S". 
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FIG. 4. (a) Effective interface potential A+ of the stripe part (dotted line) and A_ of the 
embedding substrate (dashed line) for the thermodynamic state denoted as "a" in Figs. |2| and ^ 
(T* = 1.1, A^* = 0.005). I, II, and III indicate the local minima of A±. The global minimum of 
A± corresponds to the equilibrium film thickness / on a homogeneous "it" substrate. A-|-(/ — > oo) 
increases linearly with the slope {pi — pg)Afi. (b) Equilibrium liquid-vapor interface profiles for 
the same thermodynamic state as in (a) and for the stripe width a ~ 27.5a for which at this value 
of T and Afi the morphological phase transition takes place (see Fig. Therefore the value of 
the functional is the same for both interface profiles shown here. The dotted (dashed) lines 
indicate the positions of the extrema of (A_). 
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FIG. 5. (a) Same as in Fig. §(a), but for T* = 1.1 and A^* = 0.025, i.e., for the thermodynamic 
state denoted as "b" in Figs. |2| and ^. The equihbrium interface profiles shown in (b) correspond 
to a ~ 6.4(T, for which at the given thermodynamic state the morphological transition takes place 
(see Fig. |3|). In contrast to the situation shown in Fig. Qthe repelled solution is bent towards the 
substrate and there is only one local minimum (I) of A+ . 
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FIG. 6. (a) Same as in Fig. |(a), but for T* = 1.155 and A/i* 0.0232, i.e., for the thermo- 
dynamic state denoted as "c" in Figs. |2| and |3| which is located on the hne "ii" . The equihbrium 
interface profiles shown in (b) correspond to a ~ 26a, for at which the given thermodynamic state 
the morphological transition takes place (see Fig. The positions of the minima II of A-|- and 
III of A_ coincide so that the repelled interface shown in (b) is flat: Irepi^) = I- where /_ is the 
position of both minima II and III. 
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FIG. 7. (a) Same as in Fig. §(a), but for T* = 1.15 and A/i* = 0.005, i.e., for the thermody- 
namic state denoted as "d" in Figs. ^ and ^. (b) displays the interface profiles for this state and 
for a = 50a. In this case the minimum III of A+ is deeper than the local minimum I because of 
T > Tpre- Therefore for all temperatures within the range rpj.e(A/i) < T < Tm^Afi), such as the 
one chosen here, the repelled solution (full line) is the equilibrium profile (see Fig. |3|) whereas the 
bound solution (dashed-dotted line) is metastable. The line Tm^Afj,) is shown as the dotted line 
"iii" in Fig. |. 
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